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ABSTRACT 


The modeling of physical phenomena is a useful tool 
for scientists and engineers in representing the important 
characteristics of the actual mechanism. A capillary pressure 
function for porous media (commonly known as the J-function) 
was given by Leverett and has been widely used in Petroleum 
Engineering. However, his model is known to correlate 
satisfactorily only the data from unconsolidated sands and 
from the same formation. It is the purpose of this study 
to formulate a new drainage capillary pressure function in 
an attempt to eliminate this restriction. 

The new model suggested here is a function of 
Saturation and involves three parameters which have to 
be obtained experimentally. Results obtained from the 
centrifuge method together with literature data were used 
to test this model. 

The model satisfactorily correlates literature data 
regardless of the wetting-non-wetting fluid combinations. 

A capillary pressure curve generated from centrifuge data, 
using this model, agrees well with that obtained using a 
modified Hassler method. These results indicate that the 


model is satisfactory and that it may be regarded as 


universal. 
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CHAPTER I 
INTRODUCTION 


In petroleum reservoirs two fluids exist, and often 
three fluid phases are present. The co-existence of two or 
more phases in a porous medium gives rise to capillary forces. 

Although the absolute magnitude of the capillary 
pressure in most petroleum reservoirs is not large, the 
effects of capillary forces are extremely important in 
petroleum production. Firstly, the original distribution 
of fluids in the reservoir rock is controlled by 
gravitationai and capillary forces. (34,35) Secondly, during 
the production process, the relative freedom of movement and 
the distribution of fluids are primarily influenced by these 
forces. And finally, capillary forces are responsible for 
trapping a large portion of the non-wetting fluid within the 
interstices of the rock. For these reasons the understanding 
of capillary forces is essential in petroleum production. 

It is of interest, therefore, to develop an equation 
which describes capillary effects in porous media. It is 
also of interest to express such an equation in a universal 


form so that any capillary pressure curve could be reduced to 


a Single curve. 

Leverett (28) gave the essential concepts of 
the capillary pressure in a porous medium. He defined the 
term “capillary pressure" and proposed a universal Capillary 
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2 
pressure function, commonly called the J-function, which cor- 
relates the physical properties of a porous medium and fluids 
and the capillary pressure. It has been noted, however, that 
the J-function is not truly satisfactory when data are 
taken from rock samples with various lithologic properties 
rather than from samples of similar lithologic type. A 
possible explanation to this problem may be that the 
J~function does not consider the wetting properties of the 
medium. 

It is the purpose of this study to devise a universal 
function for drainage capillary pressure-fluid saturation 
relationship. One of the parameters involved in the proposed 
capillary pressure model incorporates the effects of inter- 
facial tension, wettability and pore-size distribution. It 
is thus possible to reduce capillary pressure curves 
taken from various rock samples to a single curve. 

A technique for using the model with data obtained 
in a centrifuge is presented. The results obtained using 
the proposed model agree well with the results obtained 


using Hassler's (19) technique. 
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CHAPTER ITI 


THEORY AND LITERATURE REVIEW 


A. Surface Free Energy, Surface Tension 
and Interfacial Tension 

At the boundary of various phases in the reservoir 
exist the forces resulting from the molecular action between 
the phases. The physical forces of molecular attraction, 
generally known as van der Waals forces, are the fundamental 
basis for explaining various boundary phenomena. The van der 
Waals forces between two molecules are proportional to the 
product of their masses and inversely proportional to the 
square of the distance between them. 

A molecule in the bulk fluid is surrounded completely 
by other molecules; thus it is attracted evenly to all direc- 
tions. However, a molecule at the boundary of any two phases, 
i.e.. gas and liquid, liquid and liquid, or-.liquid and solid, is 
surrounded by the molecules of the other phase also, thus the 
net forces acting on the molecule result in unbalanced 
attraction. As a result, the boundary surface tends to 
contract to a minimum possible area. This spontaneous 
contraction of the surface indicates that free energy is 
involved. 


In order to create a unit area of new boundary, 
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4 
work comparable to the number of molecules in the unit area 


is required. It may also be said that a certain amount of 
work is required to bring a molecule from the bulk of a fluid 
into the interface against the unbalanced molecular forces. 
The amount of work necessary, in ergs, to move molecules into 
the boundary to create a unit area of new surface is referred 
to as "the surface free energy" of the substance. (19) 

The term "surface tension" is more frequently used as 
the boundary surface tends to contract and behaves as if it 
were in a state of tension. However, the surface tension 
is really a measure of the surface free energy. It is defined 
as the force in dynes acting along the boundary surface of the 
fluid perpendicular to any 1 cm in length. Thus the surface 
tension is equal, both numerically and dimensionally, to the 


surface free energy. When the two phases are both liquid, 


the term “interfacial tension" is used. 


B. Wetting 


Where two immiscible liquids or a liquid and a solid 
are in phase contact, the surface molecules of each substance 
are attracted across the interface by van der Waals forces. 
This attraction which is called the "energy of adhesion" or 
"adhesion tension", differs from interfacial tension, which 
is a property of the interface and is the differential 


resultant pull away from it. (29) 
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For an oil-water-rock system, such as a reservoir, 


adhesion tension is defined as follows: 


or 
Y ey 
cos 60 = nae es (II - 2) 
Ywo 


The adhesion tension, or the magnitude of the contact 
angle, is the measure of the wetting properties of.the solid 
by the fluids. Equation II-2 states that if adhesion 
tension is zero, that is to say, if 9 is 90°, the inter- 
facial tension Yso = Ysw, This indicates that the oil and 
the water have equal affinity to the solid surface. If, 
however, YSO > Ysw, cos © would be positive and the contact 
angle would have a value between 90° and zero; then water 
would preferentially wet the solid surface. Conversely, if 
Yso < Ysw, © would have the angle between 90° and 180°; thus 
oil would wet the solid surface. In short, the magnitude of 
the interfacial contact angle serves as a definite criterion 
for the wetting of the solid by the fluids involved. ‘4) 

Contact angle of a solid has been measured directly 
by various workers. (4723/50) A smooth surface is necessary 


in this measurement. Most solid surfaces, however, do not have 
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(f = 12) 
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perfect smoothness. Wenzel (52) introduced a roughness 


factor to correlate the actual contact angle to the observed 


angle: 

cos ©' = r cos 0 (II - 3) 
where 

oO" > measured contact angle 

<) ; actual contact angle 

r : roughness factor = actual surface/ 


geometric surface 

However, the direct measurement of contact angle 
of a porous medium is not practical as the surface is normally 
extremely rough due to the occurance of pores. Various 
authors have reported that contact angle of a porous medium 
can be calculated from the displacement (1,3,45) and the 
imbibition (5-17) properties of it. 

The importance of rock wettability and interfacial 
tension in oil production has been discussed by several 
authors. (14/37) at has been found that the oil recovery 
from preferentially water-wet rock is significantly greater 
than those from preferentially oil-wet rock. (24736) rt has 
also been found that the wetting property of the reservoir rock 
is primarily determined by the natural surface-active substances 


in the reservoir fluids. (4) Thus, greater oil production 
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efficiency is expected by adding synthetic surface-active 
substances to control the reservoir wettability during a 


water flood. (44,50) 


C2 Capillary Pressure 


The fundamental equation for the capillary pressure 


in a single straight capillary tube is given by 


2Y¥wo cose 


ee WOOL  e a 
EC dey. ogh (II - 4) 
ie 
where 
Phat Po oe Pe Capillary spressure, dynes/cm2 
AP > Py = Por gm/cc 
© : wYvadius of the capillary tube, cm 
h «= height of the water column, cm 
Equation II - 4 states that the capillary pressure in 


this system is proportional to the adhesion tension and is 
inversely proprtional to the radius of the tube. It also 
states that the capillary pressure is proportional to the 
height of the column of water. The analogy with a bundle of 
capillaries has often been used in the study of a porous 
medium as the reservoir fluids exist in spaces which are 
Capliidany insize.. 

Leverett's (28) proposed J-function is a dimensionless 
function correlating the physical properties of the porous 


medium and the fluids with the capillary pressure in an 
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attempt to express capillary pressure data obtained from 


different samples in a generalized form: 


3 (Sq) = aan [= (II - 5) 


Some investigators include contact angle in the J-function. (2) 


ee eee ey (II - 6) 
Ycos® 2 


Brown‘?). has shown that the J-function gives a good correlation 
for the capillary pressure data obtained from unconsolidated 
sands and for the data obtained from the same formation. 
However, a poor correlation is attained for data obtained 
from different formations. Thus, he has suggested that the 
use of the J-function be restricted to specific lithologic 
types from the same formation. 

A general expression for capillary pressure as a 


function of interfacial tension and curvature of the inter- 


face is due to Plgteaus!” 
Po=y (+4) ee 7) 
2A at 2 


where ry and ro are principal radii of the interface 


curvature. It is apparent from Equation II - 7 that the 


smaller the radii of the interface curvature, the greater the 
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9 
capillary pressure. Since the magnitude of the radii depends 
on the texture of pores and the saturation of the fluids 
present, it can be said that the specific pressure difference 
across the interface may be given by a certain relative 
saturation of the fluids present in the space. That is to 
say, the capillary pressure can be expressed as a function of 


Saturation, 


i, Saturation History 


It has been known that the capillary pressure for 
a porous medium is not single-valued. There is a significant 
difference between the drainage and the imbibition capillary 
pressure for the same system. McCarde11 ‘32) has explained 
this phenomenon using a capillary system whose diameter 
changes from small to large to small. From this analysis 
he has shown that the wetting fluid saturation from the drainage 
process is considerably higher than that from the imbibition 
process when both capillary pressures are the same. It is 
thus concluded that the capillary pressure-saturation 
relationship is dependent on (1) the textural properties 
of the pores, (2) physical properties of the fluids and the 
solid, i.e. wettability of the solid and the interfacial 
tension between the fluids, and (3) the history of the 


saturation process. 
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E. Methods of Measuring Capillar 
Pressure on a Small Core Sample 


Leverett (28) was the first investigator who 
conducted a capillary pressure experiment on a porous medium. 
His experiment was done on an unconsolidated sand pack by 
"height-saturation method". The principle of this method is 
to let wetting and non-wetting phases come to equilibrium 
in a vertical sand column and measure resultant saturation 
at a number of heights. The experiment was done in such a 
Manner that drainage and imbibition capillary pressure 
curves were defined. 

The height-saturation method, however, requires a very 
long sample. Since naturally-occurring rock samples available 


for laboratory experiments are usually small, this method is 


not appropriate for such samples. For small core samples, 
four basic methods have been devised. These methods are 

(1) membrane method, (2) the mercury injection method, (3) the 
centrifuge method, and (4) the dynamic method. These methods 


are described and discussed in this section. 


1. Membrane Method 

Welge (51) proposed the use of a semi-permeable 
membrane for measuring the capillary pressure. A membrane 
is selected with uniformly distributed pores of such size that 
the displacement pressure of the membrane for the selected 


displacing fluid is greater than some predetermined maximum 
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pressure of investigation. The membrane is then saturated 
with the displaced fluid. Then some porous material (tissue 
paper, fine powder, etc.) saturated with the displaced fluid 
is placed on the membrane. This porous material ensures the 
phase continuity between the displaced fluid in the core and 
that in the membrane. The sample is then placed on top of 
the porous material. The displacing fluid is confined above 
the membrane and surrounds the sample. The pressure is 
applied to the displacing fluid in a stepwise fashion with 
small increments, and at every pressure the system is allowed 
to approach a static capillary equilibrium. As the pressure 
increases, a portion of the displaced fluid is expelled from 
the sample through the membrane and the volume displaced is 
recorded to determine the saturation. Thus, a relationship 
between saturation and the capillary pressure can be obtained. 


(49) 


Fritted glass disks{9,r12,33) porcelain, ceolle= 


phane, (43) 


among others, have been used successfully as the 
semi-permeable membrane. The maximum pressure that can be 
applied to this system is the displacement pressure of the 
membrane. Due to the limitation of the membrane material 
available, the maximum pressure is not more than 100 psi so 
(7) 


far as is known. 


A different method has been devised by Hassler and 


Brunner (19) using the same principle. In their technique, 


vacuum is used rather than pressure. The wetting fluid is 
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12 
withdrawn from the sample through a membrane by vacuum. 

Among existing methods, the membrane method is considered 
to reproduce most accurately the capillary displacement in the 
reservoir. This is so because of the slowness of the approach 
to capillary equilibrium and because any combination of fluids 
may be used with the membrane method. The latter point is 
especially a marked advantage as it has been recognized that 
the use of reservoir-like fluids is preferred because the 
properties of pore space are related to the physico-chemical 
properties of reservoir fluids. This method, therefore, is 
considered to be the standard technique against which results 
from other methods are compared. The disadvantage of this 
method is the relatively slow approach to static capillary 
equilibrium after pressure has been applied. '2) Complete 


determination of a capillary pressure curve may take weeks. 


an Mercury Injection Method 


The mercury injection method was originally devised 
as a pressure porosimeter by Drake and Ritter(41,42) 45 
measure the pore size distribution of catalysts, and was 
modified by Purcel1 ‘49) for the determination of the capillary 
pressure of porous material. Purcell's equipment essentially 
consists of a fixed-volume sample chamber and a mercury 


displacement pump. 


The principle of this method is to inject mercury into 
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18 
the sample under a number of successively greater gas pressures 
applied above the mercury. As the sample chamber has a fixed 
volume, the amount of mercury injected into the sample at each 
pressure applied is equivalent to the volume of mercury 
added to the system by means of the mercury displacement pump. 
Thus, with the knowledge of the total pore volume of the 
sample, the mercury saturation can be calculated, and a capil- 
lary pressure curve is obtained. It is necessary, however, to 
correct the results as the cell expands and mercury compresses 
at high pressure. 

In this system, the mixture of mercury vapor and 
residual air is the wetting phase, and the mercury the 
non-wetting phase. In order to correlate the results from 
this method to those from the membrane method, Purcell (49) 
has suggested a constant conversion factor independent of the 
character of the reservoir. The factor is calculated by 
analogy to capillary tubes and by assuming constant surface 


tensions and contact angles of mercury and water. The factor 


is 


=e - — = = 5 (II - 8) 


where 
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14 
Pong, Pow > Mercury and water-air capillary 


pressures, respectively. 


YHg > surface tension of mercury,= 480 dynes/cm 
Yw : surface tension of water,= 70 dynes/cm 
84g + Contact. angle of mercury,= 140° 

Ow : contact angle of water,= 0° 


The contact angle and the surface tension of mercury used in 
Equation II - 8 were obtained by Drake and Ritter (42) by 
averaging the values obtained on a large variety of solid 
materials. 

There is, however, some doubt as to the validity of 
Purcell's assumptions. Henderson(21) has shown that the 
surface tension of mercury varies with the curvature of the 
interface. Pickell, et aL, have reported that mercury is 
not strongly non-wetting; rather, it partially wets the rock 
surface. These two reports suggest that Purcell's assumptions 
of constant contact angle and surface tension may not be 
valid. 

There have been some modifications as to the conver- 
sion factor suggested. One author says that the factor 
should be calculated on the assumption that the mean cur- 
vature of an interface in rock is a unique function of fluid 


saturation. '2) In this case, the factor is 
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pat = R, pe wnyals ea = 6.87 (hy =f9) 
where 

Rp ? mean curvature of the interface 


Dumore and Schols (19) have used Wenzel's roughness 
factor of 1.3 for the Bentheim sandstone-mercury-air system 
and have suggested that the apparent contact angle for this 
System can be assumed zero. This suggestion has also 
resulted in the conversion factor of 6.87. Brown ‘7) has sug- 
gested that Purcell's theoretical value for the factor, 5.25, 
is a limiting value and that a higher value should be used 
as the water surface is sensitive to contamination, aeipanys 
it quite possible that the figure 70 dynes/cm is a maximum. 
He used different conversion factors ranging from 5.4 to 8.3 


and obtained excellent results. 


Pickel1 (38) also suggests that there is a tendency 
particularly in shaly sands, for mercury residual saturations 
to be excessively high when high initial mercury saturations 
are used. This may suggest that different pore structure 
may result when mercury and its gas are substituted for the 
reservoir-like fluids, especially when a clay-like substance 
is present. Rose (40) suggests in his discussion on Purcell's 
method that hysteretic uncertainties may be encountered 
when this method is used for the study of imbibition unless 


the porous medium is inert. That is to say, in the case of 


@aeadpuoy 2 deat: ‘bape eit” tied bide bees | 
moteye tis~yruszomsenosesané miordaed ods e3 tt ae 
suf) 10d “inne toethod Snavadun sid sul beteeggue Swit Gas ; | 

oats 2a nolvesppes 22cr -oxes Bomuece oa nto man 

~pus ead "avoir .18.3 te s9tos% aolszavnos ong geile 
28.2 \20d0n2 ot 10? aylev Lsatdetosss 2! Liesmwd Fads besesp 

baap sd Stucna sufny shdgtt € ders bamustilay worsimit eek 

prism nines ot svifiesee ef spsiswe tatew dt as 

tisindenm. 6 2 mp yennyNlOk sepa? ‘elie Sede alddweog ody 42 
£.8 ut bh. 2 mer? pkipmes dsegce? aalevevnos treast2ib Boag os ~ 

ettunsy tnetisske featsdde Bas. 

yonabne: oct avedd dade esasgeue ovfs S) pteastad 
@foiserves: [evbiees yuboxam yo? .2hmee ylede at yltetnotsisq m . 
enoiteruyse yrcotom (niSing siyint nedw eid, yleviessoxe ed 03 if 
suioets atu sneaaREb fads tesppow alr sid? bese ozs 


16 


an inert porous medium the character of the interstitial 
spaces is unrelated to the physico-chemical properties of 
the’saturating® fluids". 

The advantages and disadvantages of this method may 
be summarized as follows. The advantages of this method 
are (1) the method is applicable to small, irregularly 
shaped samples, (2) a complete capillary pressure curve can 
be obtained in 30 to 60 minutes, and (3) the range of 
pressure is 5 to 10 times that of the membrane method. The 
disadvantage is the permanent loss of the sample due to the 


unrecoverable mercury. 


3. The Dynamic Method 


The dynamic method was devised by Hassler (18) 


Gr 
relative permeability studies and was modified by Brown 7) 
for the determination of dynamic capillary pressure. The 
apparatus consists of two semi-permeable membranes, one at 
each end of the sample, and completely separate flow systems 
for the two phases involved. 

Brown has used gas as the non-wetting and oil as the 
wetting fluid. The membranes are saturated with wetting 
fluid and are placed in the apparatus at each end of the 
core. The phase continuity of the oi] in the core and mem- 
branes is established by some porous material placed between 
the membranes and the core. The oil enters and leaves the 


core through the membranes, and gas through the grooves in 
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the faces of them. The central part of each membrane is 
isolated from the other part for the determination of oil 
pressure. The difference between the pressures in oil and 
gas at the inflow end is the capillary pressure. The fluids 
pass through the sample at such a rate that the pressure 
drop in the two phases across the core is equal. This is 
accomplished by adjusting the gas flow rate. Thus, 
simultaneous steady-state flow of two fluids can be estab- 
lished by regulating the quantity of each fluid entering the 
core, and a complete capillary pressure curve can be 
obtained. 

There have been several modifications of this method 
made by various investigators. Brownscombe, (8) et al., 
permitted counter flow, but experienced plugging of the 
sample. Care must be taken, therefore, to avoid possible 
plugging and the distortion of the sample in this system. 
Leas, (27) et al., eliminated the use of the outlet end disk 
and permitted gas flow while liquid saturation remained 
constant in the core. Although there still remains some 
uncertainty as to the fundamental difference between the 
static and the dynamic capillary equilibrium, Brown's 
investigation shows good agreement between the data obtained 
from this method and the membrane method for a limestone and 
a sandstone in an oil-gas system. 

The advantages and disadvantages for this method may 


be summarized as follows. The advantages are (1) capillary 
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pressure is equal throughout the core, thus the fluid 
saturation is also uniform provided the core is homogeneous, 
and (2) the end effect can be eliminated, as the capillary 
pressure is equal at both ends of the core. These two 
points result in a marked advantage for this method, 
namely jthat (3). the.dinect measurement of, isatunati.on ands the 
capillary pressure is possible, while in other methods only 
average saturation and the capillary pressure at a certain 
cross-section of the sample are obtained. The disadvantages 
are (1) relatively slow approach to capillary equilibrium as 
in the membrane method, and (2) the difference between 
dynamic and static capillary equilibrium has not been fully 


appreciated. 


4. The Centrifuge Method 


Hassler and Brunner (19) have reported the application 
of a centrifuge for measuring the capillary pressure on a 
small core sample. This method has been used by a number 
of investigators (14,26,31,46) go, capillary pressure, 
connate water saturation, and wettability studies. 

The centrifuge method makes use of a greatly 
increased drainage effect by submitting the saturated cores 
to an increased gravitational force created by the high 
normal component of centrifugal acceleration. The sample to 
be tested is rotated at a number of successively greater 


speeds, the speed being held constant until the capillary 
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equilibrium is attained at each speed chosen. The sample 
holder is designed in such a way that the amount of fluid 
displaced is collected in a pipette and can be read by the 
Operator with the aid of a stroboscope. The average satura- 
tion of the sample can thus be determined with the knowledge 
of the total pore volume. The speed of rotation is converted 
into a suitable force unit, thus a complete capillary 
pressure curve can be obtained. The capillary pressure 
curves obtained by the centrifugal method, however, are not 
a complete description of the capillary pressure versus 
saturation relationship in a porous medium, for the procedure 
does not measure directly either the saturation or the 
capillary pressure. 

The basic theory for the computation of the capillary 
pressure has been shown accurately by Hassler. (19) tf a 
cylindrical core of length L, containing a liquid of density 
0 is subjected to an acceleration g, and if the capillary 
pressure at the outer end of the core is zero, then 
the capillary pressure at height h above the lower 
end of the core will be pgh. Since the capillary pressure 
can be considered as a function of saturation, fluid satura- 
tion, conversely, can be expressed as a function of the 
capillary pressure. By definition, the average saturation 


obtained by this method is given by 
1 L 
6 = — f S(pgh) dh (II - 10) 
igi) 
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20 
If we define the capillary pressure at the inner end of the 
core (closer to rotary axis) as Z and at any point in the 


core as x, Equation II - 10 can be simplified as 


a Z 
zai= f  S(x) dx (II - 11) 
O 
where 
x = ogh 
and 
Z= ogh 


From Equation II -1ll, it follows that 


siz) = & (28) (tateees 13) 
az 


The saturation for capillary pressure Z can thus be obtained 
by graphical differentiation of the curve ZS versus Z. 

The centrifugal acceleration, however, is propor- 
tional to the distance from the centre of rotation, and will 
change along the length of the sample. Thus, the equivalent 
equation to Equation II - 11 for the centrifugal field is 


more complicated. It is given by Equation II - 13: 
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vA 
0, 
S = OS its S(x) dx 
ZS = cos 5 "f (II - 13) 
Dee aes BR | 
O Z 
where 
p 
Be ae 
iach 5 al ald ~ x,°) 
x =— w (x5? - r’) 


cos a = r/o 


The solution of equation II - 13 for S(x) involves 
very tedious successive approximations. Hassler (19) has 
suggested, however, that a constant acceleration can be used 
along the core without introducing serious error provided the 
ratio r4/Y2 is sufficiently large. If such is the case, 
Equation II - 12 can be used to determine saturation. He 
suggested the minimum value for cos ato be 0.7. Slobod, 
et al., have used the average acceleration measured at the 
centre of the core and cos% as 0.85. The capillary pressure 
in psi at height h in the core at speed n (rpm) with the 
average radius of rotation r is then given by the following 
equation: 

2 


P= 1.578 x 107/p: n* rh fread 5 
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Szabo '47) has devised a technique for measuring 
imbibition capillary pressure in a centrifuge. A sample 
having a water saturation close to the connate water 
Saturation is placed in the imbibition cell, and then a 
semi-permeable membrane permeable only to water is placed 
on. the face ofvthesinner.endhofethe-core; sAfter the 
continuity of water in the core, membrane, and the water 
reservoir located on the opposite side of the membrane is 
ensured, the cell is filled with oil. The centrifuge is 
then started at a speed less than the maximum attained 
during the displacement test. After the capillary equilib- 
rium is reached, the centrifuge is stopped to determine 
the saturation by weighing the sample. 

While imbibition takes place against gravitational 
acceleration in the membrane method, in Szabo's method, 
the centrifugal acceleration works in favour of imbibition. 
Szabo also assumes that significant re-distribution of 
fluids in the core does not occur while the centrifuge is 
stopped for the determination of saturation. 

Hoffman (22) has reported that the time required for 
determining the capillary pressure can be further reduced 
by using a constantly-accelerated centrifuge. A centrifuge 
with a speed control system for increasing centrifuge 
speed at a constant rate is used in his study. Centrifuge 
speed, average saturation, and capillary pressure are 


expressed as a function of time, and by proper rearrangement 
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3 
of Hassler's equation, saturation at the end of the core 
closer to the rotary axis can be evaluated ata given time. 
A complete capillary pressure curve can be obtained in 
approximately 6 hours by this method, while Hassler's equip- 
ment requires at least twice as much time as Hoffman's. 

The advantages for the centrifuge method are (1) the 
greatly reduced time required for obtaining a complete capil- 
lary pressure curve, and (2) high capillary pressure created 
by proper design of the equipment. The disadvantages for 
this method are that (1) the capillary pressure saturation 
relationship cannot be obtained directly, and that (2) the 
imbibition curve for pressures greater than zero and the 
drainage curve for pressures less than zero cannot be 


obtained by this technique. 
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CHAPTER III 


A CAPILLARY PRESSURE MODEL AS A 


FUNCTION OF SATURATION 


A. Capillary Pressure Model 


In order to carry out the analysis, an assumption 
was made that the capillary pressure can be regarded as the 
independent variable which determines the saturation at a 
given cross-section within a porous medium. It was also 
assumed that the porous medium is homogeneous. The func- 
tional relationship between the capillary pressure and 
saturation was such that it meets the following criteria: 

(1) The area under the capillary pressure versus satura- 
tion curve must be finite, as the area represents the 

work done in creating a new boundary surface. (28) 

(2) The capillary pressure is finite near S, = l. 

(3) The slope of the capillary pressure curve near 
Sy = Syy approaches minus NEL ys 

(4) Any functional relationship must agree with experi- 
mental evidence. 

Generally, two types of drainage capillary pressure 
curves are found among the existing data. Curve 1 in Figure 


TII - 1 is believed to represent the typical capillary 


pressure curve of a very thin porous disk, (48) while curve 2 
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Typical Capillary Pressure Curve 
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is often obtained in laboratory measurements. Note that 
there is a discrepancy between the two curves in the region 
near S, = 1. In this region, curve 2 shows inflection, 
while curve 1 does not. It is believed that for a homo- 
geneous porous medium, the capillary pressure curve would 
behave like curve 1, as the slope of the curve during 
drainage is a measure of the pore size distribution. The 
more uniform the pore size, the more nearly horizontal or the 
flatter the capillary pressure curve, (13) 

It should also be pointed out that the capillary 
pressure data obtained using the membrane method have often 
been misinterpreted. The results obtained from this method 
do not represent the true capillary pressure versus saturation 
relationship; rather, they represent the relationship 
between the capillary pressure at a particular cross-section 
of the sample versus average saturation. Thus,°serious error 
may result in determining displacement pressure and both 
Saturation and capillary pressure in the relatively low 
pressure region. 

Therefore, it seems reasonable to select curve 1 in 
Figure III - 1 as the basis for a capillary pressure model. 
It is also believed that this type of curve can model data 
from naturally occurring rock samples, for the data in the 
region near S, = 1 usually have little practical importance 


for most engineering purposes. 
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B. Construction of the Model 

It is generally accepted that the capillary pressure 
curve is discontinuous at the displacement pressure, Pg, 
below which saturation is unity. It is also generally 
accepted that there exists some minimum wetting phase satura- 
tion for a porous medium. Therefore, it seems reasonable to 
consider the domains of the capillary pressure and saturation 
as given below in order to describe the continuous portion of 


the curve mathematically: 
Pa So cae Se 


and 


The displacement pressure and the connate water saturation, 
Swi, are believed to be unique characteristics for a system 
of a given..combination of fluids and solids. (20) 

Fayer and Sheldon(16) have suggested that, for a 
water-wet medium, the slope of the capillary pressure curve 
approaches minus infinity as saturation S decreases to zero, 


where 
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This is consistent with their condition that dP./dS should 
behave like 1/S near S equal to zero. On this basis it seems 


reasonable to postulate that 


1 
SS ae 1 ee ee 
2 ( ) 
which, upon integration, yields 
Po =o ln So + PA a ae) 


(e 


The capillary pressure may be normalized by defining 
To = nc (IV tx 4) 


where o is the normalizing parameter for capillary pressure 
which incorporates the effects of interfacial tension, 
wettability, and pore-size distribution. Both Pg and o have 
to be determined experimentally. 


The model now can be written as 
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C. Application of the Proposed Model to the 
Data Obtained Using the Centrifuge 


The proposed capillary pressure model (Equation 
III - 5) has to be modified in such a way that the data 
obtained using the centrifuge (1.e., average saturation and 
rotation speed) can be used. The centrifugal acceleration 


is defined by 
ao = rw? Gongs), 


The capillary pressure in centrifugal acceleration is 
given by 


Pai= Ape (R = h/2)h (III - 8) 


Note that it is assumed that the capillary pressure at the 

outer end of the core is zero. This assumption is justified 
by the good agreement between the results obtained using the 
(19) 


membrane and the centrifuge methods. The capillary 


pressure profile within the core is shown in Figure III - 2. 


By definition, the average saturation of the core is 
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At some critical rotation speed, Wor the capillary 
pressure at the outer end of the core reaches the displace- 
ment pressure. Until: such rotation speed is reached, the 
water saturation remains unity. After the critical rotation 
speed is exceeded, the location at which Py = Pa is assumed 
to move along the core and to approach h = O. 

The displacement pressure is a unique property of the 
porous medium, and has only one value for any particular 
specimen, ‘209) Its value is a measure of the degree of rock 
wettability, the oil -water interfacial tension, and the 
diameter of the largest pore on the exterior of the rock 
sample. Therefore, if the largest pore size is uniformly 
distributed throughout the core, as it is expected to be in 
a homogeneous porous medium, the above assumption should be 
valid. Figure III - 3 shows the saturation profile in the 


core for some w>w where h* is the location in the core at 


(all 
which Po = Pg for any w>w>. 
From the derivation cited in Appendix A, the 


following expression results: 


Z - 
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Capillary Pressure Profile Within a Porous 


Medium Being Rotated in a Centrifuge 
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Saturation Profile Within a Porous Medium 


Being Rotated in a Centrifuge 
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Equation III - 11 could not be solved analytically; 
therefore, two methods were employed to solve the equation, 
eee 

(1) Numerical Integration Method 
(2) Approximation Method 

The denominator within the integral sign of Equation 

III - firanaves in a linear fashion in the integral limits. 


Thus, the midpoint may be used to approximate this term, or 
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Note that from Equation III - 8 
2. Pa 
h* = R -, /R* - 5 (III - 13) 
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The above equation was used as the Approximation Method. 
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CHAPTER LV, 
EXPERIMENTAL APPARATUS AND PROCEDURE 


In order that the capillary pressure in a porous 
medium be examined experimentally, a total of six sandstone 
cores were obtained from the Department of Chemical and 
Petroleum Engineering, the University of Alberta. These 
samples were about one inch in diameter and between one and 
one and a half inches in length. 

Prior to the capillary pressure tests, the porosity 
and the permeability of the samples were determined using a 
specially constructed Boyle's law porosimeter and a liquid 
permeameter. The porosimeter consisted of a core holder 
and an expansion chamber connected to a gas supply and a 
vacuum pump as shown schematically in Figure IV - 1. By 
calibrating with aluminum cylindrical blanks of known solid 
volume, the matrix volume of the cores was determined and 
porosity determined.. The permeameter consisted of a core 
Holder, Liquid) injection piston, flow-rate regulator, and a 
precision-bore capillary tube for flow-rate measurement. 
The permeameter is shown schematically in Figure IV - 2. The 
samples were saturated with the wetting fluid using a core 


gsaturator shown schematically in Figure IV - 3. 


A. Drainage Capillary Pressure Cell 


Two capillary pressure cells were constructed for the 
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determination of the drainage capillary pressure on a small 
core sample in the centrifuge. The centrifuge method was 
chosen primarily for the reduced time required to obtain 
data and for its reliability and good reproducibility. Ini- 
tially, the membrane method was considered for use along 
with the centrifuge method, but due to unavailability of 
membrane material with suitable pore size, the membrane 
method was abandoned and only the centrifuge cell was con- 
structed. A semi-permeable membrane was considered for use 
in the centrifuge cell, but this design was abandoned for the 
Same reason and the cell was constructed without a membrane. 
Figure IV - 4 shows the cross-section of the cell in actual 
dimension. The cell consisted of a core holder and a pipette 
for determining the volume of displaced fluid. A 3/8 inch 
diameter precision-bore (I.D. tolerence + 0.0004 inch) 
glass tube was used as the pipette. 

The core holder body was made of anodized aluminum, 
while the core support disk and the bottom cap were made 
of stainless steel. The arm and the tube body were made of 
mild steel, in order to locate properly the centre of 
weight. A stainless steel coil spring was placed in the 


core holder to stabilize the sample. 


Be The Apparatus Assembly 


An International Model U centrifuge was used in this 


work, operating at various speeds up to 2,200 rpm with a Model 
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Drainage Capillary Pressure Cell 
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40 
976 centrifuge head (Trunnion, 4 place). A window was 
provided on the lid of the centrifuge to facilitate obser- 
vation of the drainage cell with the aid of a calibrated 
stroboscope (1540-Pl-General Radio). 

After capillary equilibrium was reached at each 
centrifuge speed selected, the water volume displaced from 
the core was recorded by means of a camera (Miranda 85 mm, 

F 3.5) using high speed film (Tri-X, ASA 400). The same 
stroboscope was. used. to provide a light source for photo- 


graphy and to determine the speed of the centrifuge. 


Cc, Experimental Procedure 
ie Preparation of Sample 


The dimensions of the samples were obtained by 
taking the average of ten caliper measurements, and the bulk 
volume determined. Prior to tests, the samples were cleaned 
ina soxhlet extractor (2) for 24 hours: wach CCl ay dried at 95° 
C for 24 hours, then heated at 500°C for 10 hours to remove 
organic material from the surface. This treatment is be- 
lieved to render the samples strongly water-wet. (11) after 
the samples were cooled in a desiccator, porosity was deter~ 
mined by Boyle's law porosimeter using nitrogen. 

Dow Corning 200 Fluid (1.55 cs) was Selected as Che 
non-wetting phase for its chemical stability and purity, and 
5% brine was selected as the wetting phase for its similarity 
(49) 


to reservoir water. The densities and the viscosities 


of the wetting and non-wetting fluids and the interfacial 
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41 
tension between them were obtained at 20°C at the Standards 
Laboratory, Department of Chemical and Petroleum Engineering, 
the University of Alberta. 

In order to saturate the samples, they were placed 
in the core saturator and kept under vacuum at l0uHg for 24 
hours, then saturated with 5% brine. Liquid permeability 


was then measured using the permeameter described earlier. 


2. Capillary Pressure Determination 


After the capillary pressure cell was assembled, a 
small amount of brine mixed with uraine dye was injected 
into the pipette using a hypodermic syringe to establish the 
initial brine level. The cell was then filled with oil, 
leaving sufficient space for the sample to be tested. The 
sample, after having excessive brine blotted from its sur- 
face, was immersed in the oil and placed on the support disk 
in the cell. A stainless steel coil spring was used to 
stabilize the sample in position. Both centrifuge cells 
were prepared in this manner and the weight of both cells 
was made identical in order to minimize centrifuge 
vibration. 

The initial oil-water interface level was recorded 
with the camera, then the cells were placed in the centri- 
fuge on opposite sides of the rotation axis from each other 
and the centrifuge started. The speed was held constant at 


a number of successively greater speeds until capillary 
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equilibrium was reached. The equilibrium was normally 
obtained in 30 minutes to 3 hours depending on the sample 
and the centrifuge speed. The stroboscope was used to 
determine the speed of the centrifuge and to provide a 
source of light for photography. The uraine dye, being 
fluorescent, made the interface clearly visible and improved 
photography. Care was taken especially to maintain constant 
centrifuge speed, as the fluids in the core tend to re- 
distribute if the speed fluctuates. After the measurement 
was completed, the centrifuge was stopped and the remaining 
brine was extracted from the sample using a soxhlet 


extractor for the material balance check. 


3. Determination of Average Saturation 
and the Capillary Pressure 


Photographs were printed magnifying the pipette to 
approximately three times its size. With the aid of 
reference lines provided on the pipette, the oil-water 
interface level was determined with calipers and the volume 
of brine displaced from the sample obtained, from which the 
average water saturation could be readily calculated. 

The capillary pressure at the inner end of the sample 


was determined using the following equation: ‘48) 


$= 1.578 x 10° ap no(R = H/2)H (iy bay 


The distance from the centre of rotation to the outer face 
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43 
of the core (or to the core-support disk) in this experiment 


was measured to be 16.314 cm. 


4. Treatment of Data 

Equations III - 11 and III - 12 were employed to 
calculate the true water saturation at the inner end of the 
sample. An algorithm to perform a least-square estimation 
of non-linear parameters, developed by Marquardt, (30) was 


utilized to estimate the parameters S and Pa. Arter 


wa? 
the parameters were estimated, the capillary pressure model 
(Equation III - 3) was used to generate the capillary 


pressure curve. (See Figures V - 10 and V - 11) 
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CHAPTER VY 
RESULTS AND DISCUSSION 


A. Testing of Model 

In order to test the model, capillary pressure data 
were taken from Calhoun's (11) paper. His data were obtained 
for two silica core samples using the membrane method. The 
displacing and the displaced fluids used were air and water, 
respectively. The lengths of the samples were not given by 
Calhoun, and therefore were assumed to be 2.54 cm as this 


was the size commonly used in the membrane method. 

Two approaches were employed to test the model: 
(1) by rearranging the model to interpret the data correctly 
and (2) by assuming that the membrane method results repre- 
sent a true saturation versus capillary pressure relation- 


ship. 


ie Approach i 


As is mentioned in a previous chapter, the membrane 
method results do not represent true saturation versus 
capillary pressure relationship due to gravity. Figure 
V-1-A is a simplified illustration of the membrane cell 


for an air-water system. 


Air pressure ae is recorded as the experimental 
capillary pressure. Note that P is the capillary pressure 
g 


at h = 0. Due to gravity, the capillary pressure increases 
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with height to a maximum of Py ei 9H, (=2)- at cht Ht ‘Conises 
quently, saturation also varies along the core as is shown 
in Figure V - 1 - B. The saturation determined experimentally 
is in fact an average saturation of the core. 

It is necessary to alter the model so that average 
saturation S, can be correlated to Pg. Equations V - 1 and 
V - 2 result from the required derivation given in Appendix 


Be 


For air pressure Py > Pa 


3 (1-S,,i) Py-P Z-P 
lear Eines d 
Sw @hewaor exp (2259) exp) (= ) 
0. 9H oO oj 
(V ~ 1) 
For Pg < Pg < Pg + PwIH 
Bok Lo: 2 Z-Pq 
Sw a Swi = Like Pq = Pg aie (0) { 1 exp (s - ); 
PwIH 
(Visa) 


The initial displacement takes place when the Sap 


lary pressure at the top of the core, i.e. 2, reaches P,. 
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d 


this reason that two equations are derived, one for P 


The air pressure at this time will be P. - P.gH. It vs for 


ye ad 


and the other for Pg < Pq- 


Using Equation V - 1 and the capillary pressure data 
for Pg > Pg, the parameters were estimated applying the 
Marquardt Algorithm. These parameters were used in the 


following manner: 


(1) S versus Py generated for Py En using Equation 
V-1 


(2) Sy versus P, generated for P, < Pg using Equation 


g a 


aa 


A complete S,, versus P_ curve was obtained by combining the 


g 
results. The solid curve (curve 1) in Figures V - 2 and 
V - 3 is the generated C. versus Pg relationship for samples 
a gol and 3 = dij .respéctively. 

It is seen that the model simulates the inflection 
of curvature for the relationship near S, = 1. This inflec- 
tion is due to the fact that che is calculated rather than 


S However, the greater degree of inflection in the exper- 


SEL 
imental data indicates that it is influenced by other factors 
than just gravity, such as the pore-size distribution (i.e. 


existence of very large pores). 


ES Approach eo) 


Here, the 7 versus a results were assumed to repre- 


sent S,, versus P.. The model (Equation III - 5) could readily 
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50 
be used to estimate the parameters without any rearranging. 


Using the data points for P._ > Pq, the parameters 


g 
were estimated and the capillary pressure curve generated. 
Generated curves are given as the dotted line (curve 2) in 
Figures V - 2 and V - 3. It is to be noted that the solid 


curve in the region where Py, > Pq represents generated data 


from both approaches. 


3. Comparison of the Two Approaches 


Excellent agreement is noted between curve 1 and 2 
in. Figures V - 2 and V - 3 and between the parameter values 
estimated (Table V - 1). Parameter values are in fact the 
Same except at P,, where this difference is within linear 
confidence limits. Separation of the two curves in the re- 


gion P. < Pq is due to the neglect of gravity in the second 


g 
approach. For the same reason, the Pg estimated in the 
first approach is greater than that in the second approach. 


The gas pressure P_ at SH = 1 obtained in the first approach 


g 
is less than P, by a constant amount of p. gH. 

The agreement between the two approaces justifies 
the assumption made in the second approach. This assumption, 
however, is not valid when the sample is relatively long. 


The agreement between the data and the modeled curve indicates 


that the lengths of the samples assumed are not unreasonable. 
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52 
4. Dimensionless Caps itary Pressure 
Curve and Conclusion 

Calhoun obtained capillary pressure data on tetra- 
decane-water and oleic-acid-water systems using the same 
samples. In these systems, the gravity effect is smaller 
than in the air-water combination because smaller density 
differences exist between the fluids. Using the second 
approach, parameters were obtained and then normalized capil- 
lary pressure (Equation III - 4) and saturation (Equation 


III - 1) were calculated using the data points for P_ > Fal 


g 

The results shown in Figure V - 4 indicate that the 
suggested model can be used effectively regardless of fluid 
combination. The parameter o seems tO incorporate the effect 
of interfacial tension, wettability, and the pore-size 


distribution. With these results, it is concluded that the 


model may be regarded as universal. 


Bs Results 


1. Experimental Results 


The physical properties of the samples are given in 
Table V - 2. The samples are all sandstone and are approxi- 
mately one inch in diameter and one and a half inches in 
length. The porosities range between 15 and 24 percent while 
permeabilities range between 30 and 430 millidarcies. The 
properties of the fluids used are given in Tableev > 3. 


The experimental data ay and Z for each sample are 
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tabulated in Tables C - 1 through C - 6 in Appendix C and 

are plotted in Figures V - 5 through Vv - 10. Three duplicate 
experimental runs were made on samples A, B, and C, and two 
runs were made on the other samples. 

During the experiments, no shifting of the centre of 
rotation or irregular vibration was observed in the centri- 
fuge. The temperature in the centrifuge did not exceed the 
room temperature by more than 5°C, therefore the properties 
of the fluids were assumed to remain constant. 

The maximum time required for the samples to come to 
equilibrium was found to be approximately 2 to 3 hours at 
2,000 rpm. In the six cores studied, no change in fluid 
Saturation was observed after three hours of centrifuging. 
This was due to the high permeabilities of the samples. 

The capillary pressure is 0 psi at the outer end of 
the core (i.e. where 100 per cent water saturation is as- 
sumed to exist) and it increases with the height to a maximum 
at the inner end of the core (h = H). The maximum capillary 
pressure at the inner end of the core (one and a half inches 
in length) is approximately 6.5 psi. The typical deviation 
in the average water saturation obtained was 3 percent 
while the deviation for pressure was 0.2 to 0.3 psi. 

Table V - 4 gives a summary of the water volume 
determinations for the cores. The total volume of water 
expelled from each core at the end of centrifuging and the 


water volume extracted from each after centrifuging are 
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compared with the pore volume. The error ranges between 0.02 
and 0.09 cc. This error is equivalent to 0.5 to 2.3 per- 


cent saturation for a typical core. 


2. Numerical Results 
(1) Modeled Capillary Pressure Curve 

The Marquardt Algorithm was employed to obtain S., 
versus P, curve from the experimental data. Two methods 
were used to apply the model as was mentioned in Chapter 
Plies Cr 

(1) Numerical Integration Method 
The Marquardt Algorithm was used in conjunction 
with a Gaussian Quadrature 40-point routine to inte- 

grate the second term of Equation III - ll. 

(2) Approximation Method 
The Marquardt Algorithm was applied directly to 

Equation, GLl = 12. 

The parameters of the capillary pressure model were 
estimated with a convergence criterion of ¢ = 107°. Both 
Numerical Integration and Approximation Methods were used for 
samples A, B, and C while only the Approximation Method was 
used for the other samples. The parameter values obtained, 
their linear confidence limits, and the standard error of 
estimate are summarized in Table V - 5 and the smoothed 
values for 5, are given in Tables D - 1 through D - 6 in 


Appendix D. Generated capillary pressure curves are shown 
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66 
graphically in Figures Vv - 11 and V - 12, where the single 
curves seén are sufficient to describe both methods. 

(ii) A Modified Hassler Method 

The capillary pressure curve for sample D, run 2, was 
obtained by applying a modified Hassler method, for which the 
Spline Fit Technique (25) was used. In this technique, data 
points are arbitarily divided into several intervals, and the 
points in each interval are smoothed with a cubic polynomial, 
making the first two derivatives of the curve on the right- 
hand side of a boundary point equal to the first two deriva- 
tives of the curve on the left-hand side of the boundary 
point evaluated at the boundary point. The differentiation 
of the resulting cubic polynomial gives d(2S_) /dz, thus 
avoiding Hassler's tedious graphical differentiation. 

The results obtained using the modified Hassler 
method is given in Table D - 7 in Appendix D. The solid 
line in Figure V - 13 represents the results obtained by the 
Approximation Method while the circles represent the satura- 
tion calculated using the modified Hassler scheme. Essen- 


tially the same results were obtained using both techniques. 


C. Discussion 


lke Experimental Results 


The first runs on samples A and B showed a large 
deviation from the other two runs. There is some basis for 
believing that this behavior may have been caused by the 


unstable wetting properties of the samples at the time of the 
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70 
first run.(6) It is expected that a number of repeated runs 
may have eventually stabilized the wettability. It is for 
this reason that the results for the initial runs for all 
the samples were not taken into consideration for estimating 
parameters of the capillary pressure model except for sample 
C whose results for the first run did not show such behavior. 

The error involved in the determination of the water 
volume in the pipette is estimated to be + 0.013 cc when the 
pipette is magnified 3 times the actual size on photographic 
paper. The uncertainty in saturation would then range 
between 0.8 and 2.6 percent. The typical deviation in the 
average saturation (three percent) and the actual error in 
water volume determination (Table V - 4) indicates that the 
average saturation was determined with satisfactory accuracy. 

The stroboscope was calibrated to + 1% tolerance. 

At 2,000 rpm, this is equivalent to 0.13 psi for a typical 
Sample. Considering the possibilities of other error sources, 
such as very slight centrifuge speed fluctuation, the typical 
deviation of pressure (0.2 to 0.3 psi), may be considered to 
be small. 

Table V = 6 shows the estimated connate and resultant 
water saturations for each sample. The maximum rotation 
speed for the centrifuge was 2,200 rpm. This was insufficient 
to drive the water saturation to the connate water saturation. 
Thus the resultant water saturations were considerably 


greater than the connate water saturations. 
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2. Numerical Results 
(i) Suggested Method Using the Model 

Relatively large linear confidence limits are seen 
for o in Table V - 5 as compared with Syj and Pay andicating 
that o is sensitive to small changes in capillary pressure 
data. The phenomenon can possibly be explained by noting 
that the samples used are not homogeneous. Though the 
standard errors of estimates ranged between 0.368 x 1072 ana 
Lemon ox 1072, they are considered to be satisfactorily small. 

The differences between the parameters obtained in 
the two methods are much smaller than the linear confidence 
limits. Good agreement between the two methods leads to the 
suggestion that due to its simplicity only the Approximation 
Method need be used. 
(ii) Modified Hassler Method and Suggested Method 

The modified Hassler method has given essentially 
the same results as the Approximation Method (See Figure 
V - 13). Use of numerical differentiation rather than 
graphical differentiation eases the application of Hassler's 


(19) considerably. 


theory 
A problem with the Spline Fit Technique is the 

selection of boundary points. As the cubic polynomial for 

each interval is dependent on the choice of boundary points, 


changing boundary points normally results in some 


change in derivatives also. Apparently, no theoretical 
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studies have been made as to the optimum choice of boundary 
points associated with this technique. However, this should 
not be too serious a problem as long as data points are 


reasonably smooth. 


3. The End Effect 

It has long been noted that where there is a discon- 
tinuity in the capillary properties at the outflow end of the 
core sample, capillarity in the core tends to draw the water 
into the core from the void. The water thus accumulates in 
the boundary, and the increased water saturation causes a 


1, (28) This behavior is 


decrease in the permeability to oi 
called the "end effect" or the "boundary effect". 

In the reservoir or large-scale laboratory experi- 
ments, the distorted zone caused by the end effect has a 
relatively small effect on the interpretation of results. 
However, in a small-scale laboratory experiment, it is very 
important to account adequately for this effect. Thus, in 
the membrane method, it is customary to place some porous 
material saturated with the wetting-fluid between the sample 
and the semi-permeable membrane to ensure phase continuity 
and thus eliminate the end effect. In the centrifuge method, 
however, the end effect has often been ignored despite the 
absence of an end-butt in most experimental set-ups. The 
question arises, therefore, if the end effect would violate 


seriously the assumption made by Hassler that Po, = O at h= 0. 
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It is clear from Leverett's explanation of the end 
effect (28) that the pressure in the water must be consider- 
ably higher than that in the oil just outside the core, in 
order that displacement take place. Figures V - 14-A to 
V - 14-C show the pressure distribution in the oil and water 
within the core in three different situations. Figure 
V - 14-A is the distribution on which Hassler's assumption 
is based. Figure V - 14-B shows the pressure distribution 
expected when the capillary displacement is taking place, 
and Figure V - 14-C shows the distribution when the system 
is in capillary equilibrium. As the oil pressure is greater 
than the water pressure at the inner end of the core, the 
Capillary pressure inversion must take place somewhere in 
the core. It is expected that the location of inversion (at 
which Pc = 0) would move toward the outer end of the core as 
the system approaches capillary equilibrium. The same effect 
is expected as the centrifugal acceleration increases. 

Strictly speaking, however, the location at which 
rae 0 would never reach the outer end of the core. Thus the 
end effect violates Hassler's assumption and may result in 
higher displacement and capillary pressures until the 
centrifugal force increases sufficiently to reduce this 
adverse effect to a negligible level. Furthermore, negative 
capillary pressure at the outer end of the core indicates 
the saturation at the location just outside this end face is 


not 100 percent oil. This would result in higher observed 
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76 
average saturation. The end effect, however, does not seem 
to affect Hassler's assumption too seriously, as the more 
acceleration increases, the more the assumption seems likely 
to be reasonable. This is further confirmed by the agreement 
between the membrane and the centrifuge method results 
reported by Hassler and the centrifuge test done by Marx (31) 
with and without an end-butt. 

At this point, it should be noted that the displace- 
ment pressure estimated using the proposed method is consider- 
ably lower than the value obtained by extrapolating data 
PIOESs Flgure V > 5 tomy 10) to Se =" 4 Minrorder tc 
interpret the results adequately the discrepancies are 
tabulated in Table V - 7 and are plotted against perme- 
ability in Figure V - 15. It seems from Figure V - 15 that a 
linear correlation may exist between the discrepancy and 
permeability; the discrepancy increases with permeability. 
Leverett (28) explains this tendency by saying that the 
greater the permeability, the greater the distance over 
which the appreciable distortion takes place. The anomalous 
result for sample E can be explained by noting that it is 
much shorter than other samples. 

As to the discrepancies found between the experimen~ 
tally obtained Pg and the estimated Pq, the following 
explanation can be made. The Marquardt Algorithm estimates 
parameters using all data points equally weighted. Therefore, 


even though the data are slightly in error in the relatively 
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low pressure region, the model would estimate a Pg which 
would likely be closer to the true Pq than that obtained by 
extrapolation. As the Pg obtained by extrapolation is 
expected to be higher than the true Pa, the lower bh, 
estimated using the model may thus be justified. It should 
be noted however, that the deviation of the estimated ea 
from that observed is not necessarily a quantitative represen- 


tation of the end effect, but could also be a result of the 


heterogeneity of the sample. 
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CHAPTER VI 


CONCLUSIONS AND RECOMMENDATIONS 


ir Conclusions 


As a result of this investigation, the following 


conclusions can be made: 


(1) 


(2) 


(3) 


(4) 


A drainage capillary pressure model has been 
suggested in this study. The model describes the 
capillary pressure-saturation relationship ina 
porous medium satisfactorily. 

The model has been experimentally verified using 
data obtained from relatively homogeneous, water-wet, 
sandstone cores, and using several different fluid 
combinations. This limited amount of experimental 
data suggests that the correlation may be universal. 
The model can readily be applied to the data 
obtained using a centrifuge. Numerical Integration 
and Approximation Methods were used. Due to the 
good agreement, the Approximation Method is recom- 
mended for practical use because of its simplicity. 
Good agreement was found between the Approximation 
Method and the conventional (Hassler) method. 

The application of the model to the centrifuge data 


has the following advantages over the conventional 


method: 
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81 
(i) Data can be smoothed, the connate water 
Saturation and the displacement pressure 
estimated, and the capillary pressure curve 
generated readily and efficiently. 
(ii) Tedious graphical or numerical differentiation 
at the data points is unnecessary. 

(6) The model estimates a lower displacement pressure 
than that observed. The estimated displacement 
pressure is probably closer to the true value. 

(7) Data obtained using the membrane method (Pg versus Se) 
Can be regarded as P. versus S,, provided the core is 


sufficiently short. 


Recommendations 


(1) To further investigate the validity of this model: 
(i) The membrane method should be used in conjunc- 
tion with the centrifuge method applied here. 
(ii) Homogeneous samples should be studied. 
(iii) The end effect should be further investigated. 
(2) Regarding this experiment: 
(i) To further reduce the systematic error inherent, 
a refrigerated type of centrifuge is desirable. 
The importance of the accuracy of the speed 
control can not be over-emphasized. 
(ii) Photographic determination of the fluid volume 


displaced was efficient and accurate; however, 
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82 
the stroboscope and the centrifuge speed should 
be electronically synchronized. 

(iii) To investigate the capillary pressure of 
samples over a broader range of permeability 
and porosity would require a higher speed 
centrifuge. 

(3) Experimental runs should be repeated a sufficient 
number of times on the same sample to stabilize 


wettability. 
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Subscripts 


Ba) is) fe} jel ex 


N 


= ih) ‘tm (0) (Ee gel fel (6) 


NOMENCLATURE 


adhesion tension, ergs/cm2 
force of gravity, cm/sec2 
height of sample, cm 


height from the outer end of a core to any 
location in the core, cm 


permeability, Darcies 
length, cm 

rotation speed, rpm 
pressure, psi 


distance from the centre of rotation to the 
outer end of a core 


Saturation, fraction 


the capillary pressure at the inner end of 


a core in the centrifuge, psi 
centrifugal acceleration, cm/sec? 


interfacial tension, dynes/cm 

contact angle 

density, gm/cm3 

capillary pressure normalizing parameter, psi 
POLosity, fraction 

angular velocity, radians/sec 


normalized pressure, dimensionless 


Gapillary on critical 
displacement 

gas 

initial 

opel 

residual 

solid 

water 
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DERIVATION OF EQUATION III - 10 
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APPENDIX A 


Given: 
Po = Apw2 (R - Byh (A - 1) 
Pq 
oes it (A - 2) 
and 
h* = R - p2 . Fa (Ae) 
Apw? 


H 
. 1 C= dq 
Sw = a ia Swi + (1 - eon) exp (- e )dh 
(A - 4) 
But 
oe for Po. > a4 
Then 
h* H p 
eh Ga = pd Sa 
6) he 
(A = 5) 


Cs 


(8 - A) 


where h* is the location within a core at which ee 
From Equation A - l, 
dP 
dh = = : (A - 6) 
Vew2 
By substituting the variable h with Po and defining 
Z = bow*(R - H/2)H (A - 7) 


the following equation will result: 


con (1 = Swi)h* (1 - Swi) 
Sree eee ae acer 
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APPENDIX B 


DERIVATION OF EQUATIONS V - 1 AND V - 2 
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APPENDIX B 
FOr Ve 2) fa: 


Pa-P 
Sy = Si + (1-Syi) exp (- —2—4) (eee) 


lo} 


But, the capillary pressure within a core in the membrane 


method varies with height due to gravity. Define 


Po = Py ede fat (Bis= 22) 
and 
Z = Pg + pywgH (B= 3) 
Case 1: Py < Pa < 4 
Pa = Py + pep hee a) 


where h* is the location within a core at which Een Py 


(B = 5) 
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By substituting the variable h with P, and h* with 


(Pg - Pg) /PwIr the following equation will result: 


Sy 7 Syi + ————— [Py - Pq + o{1-exp(- ——} 
PwIH 
(Bees 7) 
Case 2: Py ZF 3 
Pe er Po — Pg 
Sw =F fe Sipe Swi) exp (- pase 3c 
aad ie 
wil Poe-Pa.aPre 
= Sui + f exp(-- d) 
H P 0 Py? 
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eS) P_-P Z-Pa 
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EXPERIMENTAL DATA 
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TABLE € = 2 =a 


EXPERIMENTAL DATA 


SAMPLE A, RUN 1 


Centrifuge Speed Capillary Pressure Average Saturation 
n Z = 
rpm psi Sw 


479 0.379 1.000 
535 0.473 1.000 
685 0.775 1.000 
762 0.959 1.000 
808 1.078 0.971 
848 Lave? 0.962 
920 1.397 0.918 
1006 1.671 0.849 
1124 2.086 0.787 
1238 27530 0.726 
1422 3.338 0.632 
1580 A, P27 0.559 
1710 4.827 0.514 
2000 6.603 0.432 
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TASLE CG —-b 32 


EXPERIMENTAL DATA 


SAMPLE A, RUN 2 


eet Speed Capillary Pressure Average Saturation 
rpm me i 
g95 0.2598 POLO e 
566 Caio. L...0'0'0 
685 On HPS e000 
835 Palo C.990 
920 AC et 0.942 
1000 aol Ors oaeL 
1065 Va B72 O8:55 
1165 2.240 Oreo, 
1245 Ys Ake ei) O57 45 
1465 3.543 0.634 
1580 4.121 0.593 
1650 4,494 0.550 
1760 Stel es O,o25 
1915 6.054 0.495 
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TABLES C = 1 =13 


EXPERIMENTAL DATA 


SAMPLE A, RUN 3 


Centrifuge Speed Capillary Pressure Average Saturation 
n Z ee 
rpm psi 
476 On 374 1.000 
648 0.693 1.000 
Sis 7.096 0.262 
890 2.208 0.952 
1000 L652 0.880 
1140 2.145 Oc 5 
1265 20042 0.726 
1420 34529 0.655 
1535 3.889 0.598 
1805 5.678 0.499 
2025 6.769 0.443 
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TABLE Co - 2 “=] 


EXPERIMENTAL DATA 


SAMPLE B, RUN 1 


Centrifuge Speed Capillary Pressure Average Saturation 
Be od w 
421 0.294 1.7000 
504 0.421 1.000 
617 0.63. 1.000 
700 On ea2 1.000 
Teeys) 08932 10:00 
824 il ay 0.976 
880 WAR 2) 0.961 
930 1,433 0.933 

1020 lel Zo 0.879 
20 20079 0.807 
1184 one 2 0.774 
1440 3.435 0.648 
1500 Bree 0.629 
1760 Bian 0.529 
2000 6.625 0.481 
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FLAT Een (lo te 2s oy 2 


EXPERIMENTAL DATA 


SAMPLE, B, RUN 2 


Centrifuge Speed Capillary Pressure Average Saturation 
ie aa 
Si2e 0.451 1.000 
620 OA6:37 1.0100 
720 O2859 1.0.00 
800 1.060 1.000 
878 ie ou) 0.995 
920 1.402 03275 
970 Laooe 0.945 
996 1.643 O27 

1044 12805 0.880 
1090 1.968 0,856 
w120 22078 0,828 
LGL2 4,304 0.614 
1746 5.049 D569 
1860 Dea 30 0.545 
2010 6.608 0.514 


2 


o00.F T£a.0 

ooo,s | 28.0 oss 
OO0.f oa0.t os 
208.0 wrt tS... are. 
aver 808.2 ose : 
ane.o gee. 1 ore. 

rse.6 eha.t ace. ae 
ona.0 208.1 TY) 
aze.0 gae.t coor. 
898.0 ato.c Ostr 


SOE <b 
eno.e 


104 


DABGE ¢C = 992) i 93 


EXPERIMENTAL DATA 


SAMPLE B, RUN 3 


Centrifuge Speed Capillary Pressure Average Saturation 
i a w 
485 0.390 1.000 
602 0.600 1.000 
700 DSi 2 1.000 
765 0.969 1.000 
840 1.169 000 
900 1.342 O3o75 
985 1.607 0.929 

1085 1950 0.867 
1180 2306 0.804 
1280 2, 7.1.4 0.746 
1360 3.064 0}. 719 
1620 4.347 0.605 
1720 4.900 0.568 
1890 & 96 0.541 
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TABLEs GC. =" 37" 1 


EXPERIMENTAL DATA 


SAMPLE C, RUN 1 


Centrifuge Speed Capillary Pressure Average Saturation 
pie ae *w 
386 0.246 1.0.00 
486 023390 1.000 
558 0.514 1.000 
616 0.626 1.000 
745 09-16 On907 
830 dh. SG Onosu 
910 lL. S66 01.835 
972 LASoe 07.825 

1022 Ushio ONBS 
1106 2e.0uke 0.732 
1240 ZIG OF..6o 
1300 2.788 0.615 
1660 4.546 0.475 
2020 6a 73h 0.430 
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TABLEMC} =. 37-22 


EXPERIMENTAL DATA 


SAMPLE C, RUN 2 


Centrifuge Speed Capillary Pressure Average Saturation 
on aay w 
ree 0.540 1.000 
665 0.730 1.000 
755 0.940 0.970 
860 e220 0.899 
965 Mere pods 0.804 

1050 Led 0.730 
1165 2.239 0.664 
1260 27029 0.616 
1485 8.000 0.544 
1700 4.767 0.485 
1870 5.768 0.442 
1980 6.467 0.421 
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EXPERIMENTAL DATA 


SAMPLE C, RUN 3 


Centrifuge Speed Capillary Pressure Average Saturation 
a sei w 
614 On622 1000 
na 5 07091 0.984 
880 AL ed OR ons] 

1000 1.650 OMA Tews 
1100 ino 6 On 2.7 
250 Za 8 0.636 
1380 cn 9 Omer 
1650 4.49] 0.474 
1830 S28 0.447 
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TABLE G- 42 = 7 


EXPERIMENTAL DATA 


SAMPLE D, RUN 1 


Centrifuge Speed Capillary Pressure Average Saturation 
ee ey, w 
400 0.262 1.000 
600 0.589 1.000 
760 0.946 L000 
910 T3356 0.960 

1000 Sr ab 0.871 
i200 2.054 0.768 
240 Hoe at 0.688 
1340 2.940 0.646 
1420 32301 0.610 
1510 Setos 0.574 
1680 4.621 moe 
1810 5.364 0.494 
1920 6.036 0.469 
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TABLE C =~ 4 = 2 


EXPERIMENTAL DATA 


SAMPLE D, RUN 2 


Centrifuge Speed Capillary Pressure Average Saturation 
Bois a w 
555 0.504 1.000 
690 0,780 1 000 
840 LAL 1.000 

1020 1.703 0.888 
165 2 222 0 .a9A 
1290 o ais 0.47 02 
1450 3.442 0.625 
1680 4.62) 0.546 
1885 5.818 0.491 
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TABDEEC = "5 3— "1. 


EXPERIMENTAL DATA 


SAMPLE E, RUN 1 


Centrifuge Speed Capillary Pressure Average Saturation 
ee a Ww 
400 05176 1.000 
600 Os8 97 1.000 
700 0.540 0.905 
830 OF759 OPTS 
960 1.016 0.666 

1080 L236 0. S586 
1200 Ley Ooty 
1350 2.009 0.467 
1450 2.47 0432 
1550 2.648 0.402 
1670 3.074 Oneg8 
L770 3.453 0.380 
1860 3.818 0.354 
L750 AO. 342 
2050 4.632 0.331 
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TABLE C= 56 = 2 


EXPERIMENTAL DATA 


SAMPLE E, RUN 2 


Centrifuge Speed Capillary Pressure Average Saturation 
site es a 
400 024.96 1.000 
600 0.397 1.000 
670 0.495 0.930 
800 0.705 0.803 
900 0.1893 0.730 

1000 1a? 0.16162 
20 Lee 0.603 
1210 1.614 02559 
oO 2.009 0.1507 
1500 2.480 0.468 
1580 Qi ion. 0.442 
1700 Chel istis 0.424 
1830 S600 0.414 
2040 4.587 0.364 
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LABLECC=-§ 6-21 


EXPERIMENTAL DATA 


SAMPLE®F, © RUN? 1 


Centrifuge Speed Capillary Pressure Average Saturation 
n 2 Sw 
rpm psi 
400 OF289 1.000 
600 OS 558 12000 
670 OG 13010:0 
800 OR957 15000 
900 Pa2t1 12000 
1000 17425 L000 
20 OS 02980 
1210 2eue9 03958 
1360 22.765 On393 
1500 3.364 OFee9 
1580 ee S74 0.794 
1700 a2 5z1 On7G5 
1830 5.007 0.739 
2040 Gisae 2 0.688 
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EXPERIMENTAL DATA 


SAMPLE F, RUN 2 


Centrifuge Speed Capillary Pressure Average Saturation 
Se ee w 
400 0.239 1.000 
600 0.538 1.000 
700 0.733 1.000 
830 13030 1.000 
960 7.8 A000 

1080 1.744 1.000 
1200 Zao 0.992 
350 2125 0.927 
1450 3.143 0.871 
1550 See OG Sako 
1670 4.169 0.803 
LO 4.684 0.774 
1860 See 0.744 
1950 5.685 0.725 
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NUMERICAL RESULTS 
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TABLE D = 1 


NUMERICAL RESULTS 


SAMPLE A 


a a ee 


Capillary Saturation Saturation (Smoothed) 

cou = Approx. Num. Int. 
Bei w ee — 

aa rene i ct ee eee ae ES 
1.096 02952 0.986 0.986 
Se 0.990 Ceo7S 0.979 
1308 Oe az 0.955 0.954 
LPS ag 0.942 0.939 07-939 
Neely 0.880 Ooo oe 0.892 
Oo e 0.9L LE 0.892 0.892 
Pore ee os 0.852 Q.852 
2.145 0.783 0.805 0.805 
2.240 Oo iSal 0.790 0.790 
PREIS 0.745 Oietag 0.742 
2.642 On 20 Doe De iisal 
S322 05655 0.650 0.650 
3.543 0.634 0.629 0.629 
3.889 0.598 O3598 0.1599 
4.121 04593 0.582 Oe oe 2 
4.494 Qeooo Of 556 Omo50 
eee ee Bs 0.525 Ome 2 Ove 1. 
eh Ao! 0.499 O02 0.509 
6.054 0.495 0.481 0.481 
Ciuc 0.443 0.458 0.458 
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TABLE D = 2 


NUMERICAL RESULTS 


SAMPLE B 


a 


Capillary Saturation Saturation (Smoothed) 
Pressure es Approx. Num. Int. 
Zo Sy Method Method 
ed Sy Sy 
Pe) A 0.995 0.987 0.986 
dae Y p 0.929 0.978 0.977 
1.402 0.975 0.968 0.968 
Weg ey 0.945 0.940 0.940 
L<607 0.929 02930 2930 
1.643 0.927 0.923 0.923 
1.805 0.880 Oseo2 0.892 
i 250 0.867 0.865 O.865 
1.968 0.856 0.861 0.862 
2.078 0.828 0.842 0.842 
2.306 0.804 0.803 0.803 
2.714 0.746 0.745 0.745 
3.063 O.. ELS 0.704 0.704 
4.304 0.614 0.606 0.606 
4.347 0.605 0.604 0.604 
4.900 0.568 0. 576 G.. S46 
5.049 0.569 0'..570 0.570 
5.730 0.545 0.545 0.545 
5.916 0.541 0.539 0.5 5.39 
6.692 0.514 0.518 0.519 
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TABLE D = 3 


NUMERICAL RESULTS 


SAMPLE C 


ae a 


Capillary Saturation Saturation (Smoothed) 
Pressure ae Approx. Num, 2ne. 
Zz Sw Method Method 
psi = ae 
Sw Sw 


0.916 0.997 0.986 0.986 
0.940 02970 0,982 0.981 
0.992 0.984 One 7 1 Ot 0 
aloo O 932 0.934 0.934 
eeu 0. 899 Gro 0.911 
teed? Oe Wey) 0.895 02595 
1.366 0.885 OC8 71 Oecy 
1206.6 0.804 0.826 0.827 
i. 559 0.825 0.820 0.821 
1.650 OTe] 0.798 0.799 
1.223 0.783 Oe Ok 0.782 
ses Re. OF 730 0.760 0.701 
T2296 Osk27 0.724 0.723 
2.018 Or 32 Orv 20 Oren 
PP aoe, 0.664 0.682 OF682 
Page Ke) 0.661 0.639 0.639 
bre Poh 0.636 0.634 0.634 
PEAS Sie 0.616 0,629 0.628 
2.788 0. 62.5 0.609 0.608 
Sarl a On G22 0.574 Dis Dales 
3.638 0.544 ESTES) 0.535 
4,491 0.474 0.489 0.489 
4.546 0.475 0.487 0.486 
4.767 0.485 0.478 0.477 
B. a4 0.447 0.452 0.452 
5.768 0.442 0.445 0.446 
6.467 0.421 0.429 0.429 
Gea. 0.430 0.424 0.424 
7.000 0.418 0.419 0.419 
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TABLE D - 4 


NUMERICAL RESULTS 


SAMPLE D 
Capillary Saturation Saturation (Smoothed) 
Pressure 
aie Approx. Method 
Z Sw 
si a 
p ce 
eons 0.888 0.890 
PP AP sys 0.794 0.789 
PRP Pao) 0.709 0:.710 
3.442 0 .\6:2'5 0.629 
4.621 0.546 0.544 
5.818 0.491 0.492 
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TABLE D = 6 


NUMERICAL RESULTS 


SAMPLE F 
ee ea er a ae en Ye ne nt ee ee 
Capillary 
Pressure Saturation Saturation (Smoothed) 
Z ae Approx. Method 
: W 
psi a 
Sw 

255 0.992 0.994 

Zeke 0.927 0. 922 

BAS ORB 0.873 

3.592 OV 815 0. 63.1) 

4.169 0.803 0.791 

4.684 OF aia 0.764 

Bare 0.744 0.744 

5,685 02725 OT 27 
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